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M.  J.  Sewell 
1.  Introduction 

The  purpose  of  these  remarks  is  to  use  elementary  mathematics  to  describe 
some  simple  connections  between  multi-valued  Legendre  transformations  and 
certain  elementary  catastrophes.  Multi-valued  Legendre  transforms 
appear  in  subjects  such  as  nonlinear  elasticity  and  non-convex  optimization. 


2.  Some  simple  Legendre  transforms 

Consider  the  general  polynomial  equation  in  x of  the  nth  degree.  This 

can  be  written,  after  normalising  the  leading  coefficient,  ar 

n . n-1  . o 

u ■ x ♦ tx  ♦ . . . ♦ bx^  + ax 

where  the  coefficients  are  then  a,  b,  ....  t,  u.  Now  express  this  in 
the  alternative  gradient  form 


dX 

U " dx  ' 

where  the  generating  potential  function  X(x)  is  defined  by 
X 


1 n*1  t n baa? 

...+■=■  xs  ♦ — xz  . 


n ♦ 1 n 3 2 

Recall  that  such  a polynomial  equation  has  either  0 or  2 or 
or  n real  roots  if  n is  even,  and  either  1 or  3 or  ... 
real  roots  if  n is  odd. 


(1) 


■ • • 
or  n 


Those  points  of  the  potential  function  X(x)  which  satisfy 


d2X 

d^7 


0 


(2) 


can  be  divided  into  two  sets,  according  to  whether  the  lowest  value  of 


r > 2 which  makes 
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(3) 


is  (a)  odd  or  Cb)  even. 

(a)  These  points  are  inflexion  points  of  X(x),  because  the  slope  of 
X(x)  attains  a stationary  local  minimum  or  maximum  as  it  passes 
through  such  a point,  by  the  Taylor  expansion. 

Cb)  These  points  are  not  inflexion  points  because  the  slope  of  X(x) 
only  becomes  stationary  at  such  a point,  and  then  resumes  its 
previous  monotonically  increasing  or  decreasing  character.  The 
simplest  example  is  the  increasing  slope  of  X(x)  = x4  on  both 
sides  of  x = 0,  where  the  lowest  value  of  r is  4. 

If  the  normalized  constant  u in  the  original  polynomial  is  now 
regarded  as  a controllable  parameter,  equation  (1)  can  be  viewed  as  a 
mapping 

u » u(x) 

of  a real  x-axis  onto  a real  u-axis.  As  x passes  through  a point 
where  (2)  holds,  the  corresponding  point  on  the  u-axis  will  reverse 

t 

direction  for  a point  of  type  (a)  (i.e.  for  an  inflexion  point  of  X(x))j 
but  for  a point  of  type  (b)  it  will  maintain  the  same  direction  which 
it  has  on  either  side  of  the  considered  point. 

The  direction  of  movement  of  the  u-point  will  also  be  maintained,  a 
fortiori,  as  x passes  through  any  point  where 
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Such  a Legendre  transformation  is  valid,  in  the  first  instance,  within 
any  x-domain  bounded  by  points  satisfying  (2),  and  there  will  be  a 
corresponding  u-domain.  Beyond  this  x-domain  boundary  (or  boundaries, 
at  either  end)  there  will  be  other  x-domains,  leading  to  different 
u-domains  and  corresponding  different  branches  of  the  dual  function 
U(u).  However,  as  x passes  through  any  boundary  point  which  is  of 
type  (b),  i.e.  one  which  does  not  specify  an  inflexion  point  of  XCx) 
and  for  which  the  corresponding  u-point  maintains  its  direction,  the 
two  adjacent  branches  of  U(u)  can  be  matched  in  value  and  slope  from 
opposite  sides  of  their  common  u-domain  boundary.  The  second  derivatives 
become  infinite  by  (2)  and  (6),  but  the  slopes  match  because  of  (4) 
and  the  contiguity  of  the  two  x-domains. 

For  this  reason  the  convention  is  now  introduced  that  any  pair  of 


-3- 


x-domains  which  are  adjacent  at  a point  of  type  (b)  are  counted  as  two 
parts  of  the  same  enlarged  domain,  and  the  matched  pair  of  branches  of 
U(u)  are  counted  as  a single  branch  over  the  single  u-domain  formed  by 
glueing  together  the  previous  abutting  pair.  By  this  convention  a 
non-inflexion  point  which  satisfies  (2)  is  in  future  counted  as  an 
interior  point  of  an  associated  Legendre  domain. 

This  leaves  the  following  situation.  Equation  (2)  is  a polynomial 
equation  having  at  most  n - 1 real  roots.  Those  roots  which  specify 
inflexion  points  of  X ( x ) divide  the  real  x-axis  into  at  most  n 
contiguous  domains  Dx1 , D^.  ....  D (say).  As  the  x-axis  is 
traversed  from  -®  to  +“,  the  sign  of  the  curvature  of  X(x)  will 
alternate  with  each  succeeding  domain,  changing  only  at  the  inflexion 
points.  Within  each  domain  the  curvature  will  not  change  sign,  being 
either  convex  or  concave. 

If  n is  odd  there  will  be  1 or  3 or  . . . or  n such  x-domains, 
whereas  if  n is  even  there  will  be  2 or  4 or  ...  or  n of  them. 

The  sketches  in  Figs.  1-4  illustrate  those  domains  for  the  cases 

n * 1,  2,  3,  4. 

1 

As  the  x-axis  is  traversed,  the  succession  of  domains  D ^ , D •••» 
will  give  a succession  of  distinct  branches  l^Cu),  L^Cu),  •••  of  the 
composite  Legendre  dual  function 

U^tu),  ...}  = W(u)  (say), 

each  branch  being  defined  over  an  associated  u-domain  , D to 
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3.  One  or  three  curvature  domains  for 


be  determined.  A good  qualitative  picture  of  the  way  these  branches 
fit  together  can  be  built  up  on  the  basis  of  three  simple  observations: 

1.  Each  u-domain  must  be  traversed  in  the  opposite  sense  to  the  preceding 

one,  and  therefore  every  pair  of  contiguous  u-domains  must  overlap  on 
tha  whole  of  at  least  one  of  them.  This  follows  from  equation  Cl) 
and  property  (a)  above. 

2.  The  slope  of  succeeding  distinct  branches  of  WCu)  is  continuous 

at  their  join.  This  follows  from  equation  C4)  and  the  contiguity  of 
successive  x-domains.  (It  is  also  true  for  points  of  type  (b),  a 
fact  already  utilised  in  the  above  convention). 

3.  The  curvature  of  succeeding  distinct  branches  of  WCu)  alternate  in 

sign,  because  this  property  holds  for  X(x)  (see  Figs.  1-4), 
and  convexity  Cor  concavity)  is  transmitted  from  one  generating 
function  to  its  dual  in  a Legendre  transformation. 

These  observations  are  employed  in  Figs.  5 - 8 to  sketch  the  qualitative 
form  of  the  composite  Legendre  transforms  WCu)  of  the  functions  X(x) 
shown  in  Figs.  1-4.  In  these  diagrams  no  attempt  is  made  to  give  the 
quantitative  effects  of  all  the  different  possible  combinations  of  the 
parameters  a,  b,  c (for  example,  to  change  the  sign  of  a in  Fig.  2 
will  tilt  the  cusp  of  Fig.  6 downwards  instead  of  upwards.)  The  purpose 
here  is  to  convey  the  qualitative  features  which  can  be  deduced  about  the 
Juxtaposition  of  the  domains  and  branches  of  the  composite  Legendre 
transform  WCu),  on  the  basis  of  the  foregoing  three  observations. 
Succeeding  branches  of  WCu)  always  join  at  cusps. 
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Fig.  7.  One  or  three  overlapping  domains  of  dual 


Fig.  S 


W(u)  = {Ur  U2,  U } to  Fig.  3 


. Two  or  four  overlapping  domains  of  dual 


W(u)  * {U  . U_.  U . U } to  Fig.  4 


r r ~ 


Finally  it  may  be  noticed  that  if  X(x)  is  any  sufficiently  smooth 
function,  not  necessarily  a polynomial,  then  as  x passes  through  any 


point  satisfying  (2),  the  two  Legendre  transforms  U^tu),  Ui+1(u) 
defined  on  either  side  by  (5)  have  continuous  value  and  slope  at  the 
considered  point.  Locally  therefore,  if  the  considered  x-point  is 
an  inflexion  of  X(x)  (i.e.  of  type  (a)),  the  two  branches  of  the  com- 
posite transform  are  defined  over  overlapping  domains  and  meet  in  e. 
cusp;  whereas  if  it  is  not  an  inflexion  (i.e.  type  f b ) ] the  two  branches 
are  defined  over  abutting  domains  and  can  be  glued  together  to  look 
like  one  branch.  The  cusp  can  never  point  straight  upwards,  nor  can 
the  composite  transform  have  a vertical  tangent,  by  (4)  for  finite  x. 

V 

3.  Some  elementary  catastrophes 

It  can  be  immediately  recognized  that  the  composite  Legendre  dual  functions 
W(u)  portrayed  in  Figs.  6,  7 and  8 are  qualitatively  the  same  as  at 
least  certain  cross-sections  of  the  bifurcation  sets  of  the  cusp, 
swallowtail  and  butterfly  elementary  catastrophes  respectively  - see, 
e.g.  §5.3  of  Thom’s  book  (1975)  for  pictures  of  the  latter. 

V 

A complete  statement  of  the  Legendre  transformation  beginning  from  the 
above  polynomial 

X(x;  a,  b,  ....  t) 

takes  account  of  the  fact  that  the  coefficients  a,  b,  ...»  t play  the 
role  of  passive  variables.  They  therefore  enter  each  branch 
U^(u;  a,  b,  ....  t) 

of  the  composite  dual  function 
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W « W(u»  a,  b,  . . . , t) 
and  have  the  properties 


(7) 


da  3a  3b  3b 


3t  3t 


for  each  U^.  In  future  3/3x  is  also  used  in  place  of  the  d/dx  of 


Section  2.  and  3/3u  in  place  of  d/du. 


A significant  feature  of  the  dual  function  (7)  is  that  it  depends  only 
on  the  coefficients  or  'controls’  present  in  an  associated  potential 
function 

V(xj  a,  b,  ....  t,  u)  = X - ux 


n+1 


n + 1 


t n b 3 

♦ — X + . . . + -r  X3 


a o 

♦ ^ x ■ ux 


[8) 


in  terms  of  which  the  starting  point  (1)  of  the  Legendre  transformation 
can  be  written 


3V 

37  “ ° * 


(9) 


In  mechanical  contexts  where  V is  potential  energy  depending  on 
displacement  x,  with  various  possible  interpretations  for  the  control 
parameters,  (9)  defines,  in  the  space  of  x and  the  parameters,  an 
'equilibrium  surface’  as  introduced  by  Sewell  (1966,  1976).  Therefore 
(1)  also  defines  this  surface. 


Two  simple  examples  of  the  dual  function  W(u)  are  now  displayed 
explicitly. 
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Example  1 : fold  -*•  cusp. 

If  V in  (8)  is  the  universal  unfolding  for  the  fold  catastrophe,  then 
the  corresponding 


Then  (1)  generates  u * x2,  and  (2)  and  (31  hold  at  x = 0 with 

r * 3,  corresponding  to  the  (horizontal)  inflexion  in  the  cubic.  The 

origin  divides  the  x-axis  into  two  domains  : x < 0 and 

0 _ : x > 0.  The  inverses  of  u = dX/dx  in  these  domains  are 

x2 

x * -u^  and  x = u^  respectively,  both  branches  being  defined  over  the 
same  domain  0u1  * 0^  '•  u > 0 of  the  u-axis.  The  associated  composite 
Legendre  transform  is  thus 

3 3 

W - (U1 (u) , U2(u)},  U = - | u2.  U2  = | u2  . (11) 


The  two  dual  generating  functions  X and  W in  this  Legendre  transform 
are  sketched  in  Fig.  9.  This  is  a special  case  of  Figs.  2 and  6. 


Fig.  9.  Fold  potential  cusp  dual 
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Evidently  either  branch  of  the  dual  function  (11]  satisfies 


This  is  the  bifurcation  set  of  another  potential  which  is  the  cusp 
catastrophe  unfolding 

PUj  u,  W)  = - |x2uE2  ♦ 3A3W£  . (12) 

Here  A is  an  arbitrary  scale  factor  and  £ is  some  other  behaviour 
variable  identifiable  as  Ax.  This  can  be  verified  directly,  or 
inferred  from  (34)  and  (35)  below.  One  of  the  controls  in  (12)  is  the 
original  one  u,  and  the  second  control  W is  the  Legendre  transform 
generated  by  the  original  fold  potential. 

Example  2 : cusp  -*■  swallowtail. 

This  time  suppose  that  V in  (0)  is  the  universal  unfolding  for  the  cusp 
catastrophe.  The  term  of  penultimate  highest  power  is  therefore  absent 
again.  (Potentials  with  this  property  are  associated  with  singularities 
called  cuspoids.)  For  X this  implies  the  partial  unfolding 

X ■=  -J-  x4  ♦ i ax2  (13) 

4 2 

as  the  starting  point  X(x«  a)  of  the  Legendre  transformation.  There- 
fore (1)  becomes 

u ■ x3  ♦ ax  (14) 

which  defines  the  familiar  smooth  folded  equilibrium  surface  in  x,  a,  u 
space  (e.g.  Fig.  11  of  Sewell,  1976,  or  the  numerous  illustrations  by 
Zeeman,  e.g.  Isnard  and  Zeeran,  1974).  Then  the  limits  of  single-valued 
inverse  branches  of  (14)  are  defined  by  (2),  namely 

3x2  - -a  . (15) 

because  (3)  holds  with  r * 3 (except  at  the  single  point  x * 0) . 
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If  a > 0 there  is  no  such  real  limit  to  x.  i.e.  the  quartic  (13) 


i 


has  no  inflexion,  and  over  tha  whole  x-axls  there  is  a unique  real 
solution  x * x(u,  a)  to  the  cubic.  This  inverse  has  the  gradient 
x « 3X/3u  of  the  unique  Legendre  transform 

U(u»  a)  * j ux(u,  a)  - ^ a[x(u,  a)]2  (16) 

defined  over  the  whole  u-axis.  The  first  diagrams  in  each  of  Figs.  3 
and  7 give  the  qualitative  picture. 


If  a < 0,  however,  two  real  turning  points  of  the  cubic  (14)  separate 
three  contiguous  x-domains 


where  (14)  gives  solutions  x^u,  a),  x^(u,  a),  x^fu,  respectively. 

We  do  not  need  their  explicit  form  here.  Substituting  each  in  turn  into 
the  same  formula  (16)  already  derived  for  the  case  a > 0 leads  to 
three  distinct  branches 


U1  (u» 

a) 

in 

0 , : 
ul 

-«  < 

u < - 

M-M1 

U2(uj 

a) 

in 

°u2  ! 

3 

L 1.  al 

l 3aj 

^ <u<-- 

U3(ui 

a) 

in 

Du3  ! 

2a 
+ 3 

r i i 

“ 3 aj 

^ < u < •*•“> 

of  the  composite  Legendre  transform 

W(ut  a)  = {U^Iut  a),  L^Cui  a),  U^fuj  a)}  . (17) 


The  appropriate  qualitative  pictures  this  time  are  the  last  diagrams  in 
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each  of  Figs.  3 and  7,  but  symmetrized  with  respect  to  both  horizontal 
and  vertical  axes  because  of  the  absence  of  the  cubic  'bias'  term  in  (13). 

Consider  now  what  happens  when  the  passive  variable  a is  allowed  to 
vary.  The  most  significant  part  of  the  surface  (13)  is  displayed  in 
Fig.  10  by  computing  its  cross-sections  a = constant  and  x = constant. 

The  latter  are  straight  lines  with  slope  lx2.  The  former  are  quartics 
with  inflexions  whose  locus,  for  successive  a-values,  is  also  shown  on 
Fig.  10.  This  locus  has  the  parabolic  projection  (15)  on  the  x,  a 
plane,  separating  D ^ (for  each  a)  from  0^^  and  D^* 

The  induced  domain  boundaries  in  the  u,  a plane  calculated  above  lie  on 
the  familiar  cusped  curve 

4a3  ♦ 27u2  = 0 (18) 

(because  the  calculation  defining  them  is  equivalent  to  eliminating  x 
from  3V/3x  = 32V/3x2  ■ 0 for  V c lx1*  ♦ lax2  - ux,  i.e.  it  is  equivalent 
to  calculating  the  bifurcation  set  for  the  cusp  catastrophe.) 

The  composite  Legendre  transform 

W - W(ui  a)  , (19) 

is  a surface  in  three  dimensions  when  both  the  active  variable  u and 
the  passive  variable  a are  allowed  to  vary.  This  is  computed  from  (16) 
for  a > 0,  and  from  (17)  for  a < 0,  and  its  most  significant  part  is 
displayed  in  Fig.  11  via  its  successive  cross-sections  a = constant.. 

Hidden  portions  are  indicated  by  broken  lines.  To  obtain  the  best  display 
in  the  Figs.  10  and  11  it  was  necessary  to  use  rather  different  viewpoints 
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and  a-ranges.  The  Key  point  is  that  cross-sections  of  Figs.  10  and  11 
with  each  given  plane  a = constant  are  Legendre  transforms  of  each  other, 
and  the  pictures  show  how  these  transforms  vary  with  the  passive  variable  a. 


It  is  evident  that  Fig.  11  is  of  the  same  qualitative  type  as  the  bifurca- 
tion set  for  the  swallowtail  catastrophe.  As  a decreases  through  the 
origin,  the  surface  retains  its  single-valued  character  only  outside  the 
cusped  curve  (18),  while  within  it  (i.e.  over  the  domain 
27u2  < -4a3,  a < 0)  the  surface  (19)  becomes  triple-valued.  The  middle 
concave  branches  in  Figs.  10  and  11  map  into  each  other.  So  do  the 
outer  left  and  right  convex  branches  respectively,  but  with  overlap  over 
the  central  cusp-shaped  domain. 


The  local  shape  of  (19)  near  the  cusp  point  u = a = W = 0 itself  is 
hidden  in  Fig.  11.  This  is  displayed  in  Fig.  12  by  removing  the  right- 
hand  halves  u > 0 of  the  upper  two  branches  of  the  surface,  and  using  a 
slightly  different  viewpoint.  Not  only  is  the  cusp-shaped  boundary  of 
the  middle  branch  (having  projection  (18))  more  clearly  seen;  but  the 
cusped  nature  of  the  vertical  section 
W ■=  W(0>  a) 

is  also  revealed.  The  latter  is  where  the  intersection  of  the  two  upper 
branches  meets  the  a-axis  u = W = 0,  which  lies  in  the  lower  branch. 
Analytically  this  vertical  section  through  the  plane  of  symmetry  is  the 
parabolic  cusp 


W « 


- co  < a < oo 


- ob  < a £ 0 , 


(20) 
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Just  by  putting  u = 0 in  (14)  and  (16).  This  is  exactly  the  type  of 
cusp  which  appears  in  the  bifurcation  set  for  the  swallowtail  (see 
Thom,  1975,  p.65). 


In  fact  these  findings  can  be  summed  up  by  the  remark  that  if  one  begins 
with  (13)  (Fig.  10)  and  computes  its  composite  Legendre  transform  (19) 
with  passive  a,  this  surface  displayed  in  Fig.  11  is  exactly  the  same 
as  the  bifurcation  set  of  another  potential  which  is  the  swallowtail 
catastrophe  unfolding 

P(£j  a,  u,  W)  = ±C5  + §A2a£3  - X3u£2  ♦ 4A4W£  (21) 

for  another  behaviour  variable  £ = Ax,  where  A { 0 is  any  scale 
factor.  The  reason  is  that  the  two  calculations  are  equivalent  as  can 
be  seen  directly  or  from  (34)  and  (35)  below.  Two  of  the  controls  in 
(21)  are  the  original  u and  a in  (8)  (for  n = 3 with  b = 0 so 
that  (9)  =>  (14)),  while  the  third  control  W is  the  Legendre  transform 
W generated  by  the  partial  unfolding  (13)  of  the  cusp  catastrophe. 


We  investigate  the  shape  of  Fig.  11  a little  more  closely.  Every  branch 
U^uj  a)  of  (19)  has  gradients  SL^/Su  = x and  SU^/Sa  = - $x2.  Every 
branch  of  the  solution  x(u.  a)  of  the  cubic  (14)  has  the  properties 
3x/3u  ■ 1/(3x2  ♦ a)  and  3x/3a  = -x/(3x2  ♦ a).  Therefore  each  branch  of 
(19)  has  second  derivatives 


Su2  3xz  + a ' 3u3a 


3xz  ♦ a ' 


32U.  2 

i m x 

3a*  3x^  ♦ a 


(22) 


and  so  the  2*2  determinant 


I 32lJi 

I 3(u,  a) 


(23) 
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on  every  branch.  Formulae  (22)  are  directly  expressed  in  terms  of  the 
inside  or  outside  of  the  locus  of  inflexions  on  Fig.  10.  It  follows 
that  the  upper  two  branches  of  Fig.  11  are  strictly  convex  in  u and 
a separately,  but  only  weakly  Jointly  convex  in  them  (this  applies  in 
particular  where  the  upper  branches  merge  in  a > 0).  The  lower  middle 
branch  U^tuj  a > 0)  is  strictly  concave  in  both  variables  separately, 
but  only  weakly  jointly  concave  in  them. 


It  also  follows  from  (23)  that  through  every  point  of  the  surface  of 
Fig.  11.  there  is  a curve  in  the  surface  along  which  the  Gaussian  curva- 
ture is  zero.  The  local  direction  of  the  projection  of  this  curve  onto 
the  u.  a base  plane  is  given  by  the  eigenvector  associated  with  the 
zero  eigenvalue  of  (23),  namely  du/d a = x = 3Uj/3a.  In  other  words 

it  is  given  (from  (14))  by  the  solution  of 
3 


t . > 

du 

da 


du 

♦ a-; u 

da 


appropriate  to  the  considered  branch.  In  fact  each  plane  section 
x = constant  of  Fig.  10  maps  into  a plane  section  u = ax  ♦ x3  of  Fig.  11. 
This  satisfies  the  differential  equation  and  so  the  dBsired  projection  is 
a straight  line.  As  x is  altered  parametrically  these  lines  envelop 
the  cusped  domain  bo  :ndary  (18),  and  their  points  of  tangency  with  that 
envelope  correspond  to  where  the  straight  sections  x = constant  in  rig.  10 
cross  the  locus  of  inflexions  there.  The  curves  along  which  the  Gaussian 
curvature  of  Fig.  11  is  zero  are  therefore  actually  straight  lines  in 
three  dimensions: 


t ^ = 


1/x' 


W - ^x4  «=  — 
4*  2772 


(24) 
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for  each  fixed  x.  The  surface  [19]  in  Fig.  11  could  therefore  have 
been  displayed  alternatively  entirely  as  a ruled  surface  in  terms  of 
these  straight  lines,  which  transfer  from  an  upper  to  the  lower  branch  of 
the  surface  at  their  tangency  points  with  the  'edge  of  regression' 

W - - |k\  u = -2x3,  a - -3x2.  (25) 

The  latter  is  displayed  in  Fig.  11  as  the  boundary  between  the  lower  and 
upper  branches  of  the  Legendre  transform.  It  is  the  mapping  of  the  locus 
of  inflexions  of  Fig.  10. 

Fig.  13  displays  (part  of)  Fig.  11  as  a ruled  surface,  generated  from 
successive  straight  lines  (24)  by  altering  the  parameter  x.  In  order 
that  the  picture  should  not  be  self-obscuring,  we  have  selected  only  a 
limited  set  of  x-values  between  0 and  2.5^  = 1.58,  and  plotted  the 
lines  from  the  same  viewpoint  and  in  the  same  'box'  as  in  Fig.  12.  The 

lines  envelop  the  edge  of  regression  (25)  and  touch  it  at  a-values 
0,  -0.15,  -0.5,  -1,  -1.5,  -3,  -4.5,  -7.5.  Fig.  13  therefore  shows  how 
the  Legendre  transform  maps  certain  straight-line  sections  x = constant 
of  Fig.  10,  in  the  range  0 £ x i 1.58.  into  straight  line  sections  of 
Fig.  11. 
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4.  A ladder  for  the  cuspoids 


The  results  from  the  two  examples  of  53  can  be  presented  as  the  chain  of 
connections  illustrated  in  Fig.  14.  This  expresses  the  hierarchy  of 


Fig.  14  A ladder  for  the  cuspoids 

elementary  catastrophes  in  a way  different  from  any  that  I have  seen  before. 
On  the  right  in  Fig.  14  are  surfaces  entirely  in  control  spaces,  while  on 
the  left  one  behaviour  variable  also  appears  in  each  potential.  For 
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example,  from  a cusp  catastrophe  potential  (12)  we  can  pass  down  the 

1 - 

ladder  to  the  fold  potential  V * - ux,  or  (after  altering  the  notation 

to  identify  the  first  two  terms  of  (12)  with  (13))  up  the  ladder  to  a 
swallowtail  potential.  Each  of  the  bifurcation  sets  may  also  be  associ- 
ated with  many  behaviour  variables  in  the  sense  of  the  general  qualitative 
equivalence  theorem  of  catastrophe  theory. 


To  give  a general  proof  that  the  ladder  of  Fig.  14  can  be  extended  up 
any  number  steps,  it  is  convenient  to  rewrite  the  potential  (8)  as 

n 


V5’  an V = n ♦ 1 


1 ^n+1  + 


r 1 rr 
1 ~ ■ 


r=1 


r r 


(26) 


The  subscript  n is  attached  to  Vn  for  emphasis,  and  when  an  = 0 we 
may  call  V the  nth  cuspoid  potential  (cf.  Woodcock  and  Poston  1974). 


The 


fold  = -jE3  ♦ a^C  is  the  second  cuspoid  after  the  trivial  K2- 


The  bifurcation  set  of  (26)  is  the  surface  in  the  a -parameter  space 

r 

obtained  by  eliminating  £ from  3V  /3£  = 0 and  32V  /3£2  = 0,  i.e. 

n n 


from 


t"  * l v1"1  ■ »• 

r-1 

and 

n 

,n-1  r , rr-2 

nC  * / (r  - 1)a  C = 0. 

_ r 

r=2 


Eliminating  a^  from  the  first  of  these,  and  rewriting  the  resulting 
equations  in  reverse  order,  gives 
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(27) 


-a  » n^'1  ♦ l (r  - 1)a/'Z, 
r*3 


.r-2 


a.  - Cn  - 1Kn  ♦ l Cr  - 2) a £r_1. 
• _ - r 

r*3 


(28) 


Thus  the  bifurcation  set  is  found  in  principle  by  eliminating  £ between 
these  two  equations. 


On  the  other  hand,  if  we  begin  with  the  potential  V of  degree  one 
less  than  (26),  and  write  it  in  the  form 

n-1 

l 

r=1 

its  'equilibrium  equation'  3V  ,j/3x  = 0 can  be  expressed  as  the  starting 
point 


Vn  ,(xj  b .....  b)  * ^-xn  ♦ £ — b xr. 

n-1  n-1  1 n L r r 


(29) 


. 3X 

b1  * 3x 


(30) 


of  a Legendre  transformation  generated  by 

b ....  b ) = ^xn  ♦ I ~ b xr. 
n-1  2 n L r r 

r*2 

Then  (30)  is 


(31) 


. n-1  nr"'  r-1 

x * 1 b x . 

r-2 


(32) 


The  Legendre  transform  Wf-b^  b^,  ....  b2)  itself  is  obtained  by 
inserting  the  successive  branches  of  the  solutions  x of  (32)  into 
W • -b^x  - X.  By  (31)  and  (32)  this  has  the  values  of 


- - ^ • X NV- 


(33) 
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The  Legendre  transform  is  thus  obtainable  in  principle  by  eliminating  x 
between  (32)  and  (33). 


| 

I 


The  pair  of  equations  (27)  and  (28)  can  be  identified  precisely  with 
the  pair  of  equations  (32)  and  (33)  if  we  introduce  an  arbitrary  non-zero 
scale  factor  X such  that 

e - Xx  (34) 

and  then  match  the  coefficients  according  to  the  formulae 

a.  = nXnW,  a « — Xn  rb,  1 st  r S n - 1.  (35) 

i r*l  r r 


This  establishes  the  desired  general  theorem:  that  the  bifurcation  set 

obtained  in  the  n-dimensional  a -space  from  (26)  is  the  same  as  the 

r 

Legendre  transform 

W - W(-b^ i bn1,  ....  b2)  (36) 

generated  with  active  x from  (31)  in  the  n-dimensional  space  spanned  by 
W and  the  br>  provided  only  that  the  trivial  scaling  (35)  is  applied 
to  the  n axes. 


In  particular  the  result  holds  for  the  cuspoid  potentials  in  which 

a ■ 0 in  (26)  and  b „ * 0 in  (29).  Woodcock  and  Poston  (1974)  have 
n n-i 

computed  the  bifurcation  sets  for  all  the  first  seven  cuspoids  (i.e.  up 

to  n * 7 in  (26)  with  a ■ 0),  in  terms  of  various  ruled  surface 

n • 

projections  parametrized  by  £ ■ constant.  The  same  pictures  are  therefore 
Immediately  available  as  explicit  illustrations  of  the  multi-valued 
Legendre  transform  (36)  generated  from  (31). 
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5.  A mechanical  Interpretation 


The  general  theorem,  or  the  ladder  of  Fig.  14,  can  be  interpreted  as 
setting  up  a certain  sequence  of  correspondences  between  different 
conservative  systems. 


For  example,  the  two-bar  frame  analysed  by  Koiter  (196B),  or  the  spring 
model  examined  by  Budiansky  (1974),  have  a total  potential  energy  of 
the  type 


1 3 
— X J 

3 


-Ife  x2 
2°2X 


V 


where  x is  buckling  deflection,  b^  is  applied  load  and  b1  is 
imperfection.  The  bifurcation  set  of  this  model  is  the  relation  between 
imperfection  and  buckling  load 


(37) 


This  is  a parabola  in  b1 , b2-space  displaying  imperfection-sensitivity 
of  the  type  first  analysed  by  Koiter  (1945)  - i.e.  a small  (second-order) 
change  in  imperfection  permits  a large  (first  order)  decrease  in  buckling 
load. 


Stepping  first  to  the  bottom  step  of  the  ladder  (not  shown  in  Fig.  14,  but 
see  Sewell  1976,  p.  169)  (37)  is  also  the  relation  between  a weight  -^b., 

f2 

suspended  in  equilibrium  from  a Hookean  spring,  and  the  complementary 
energy  b^  of  the  spring.  For  when  the  extension  of  the  spring  is  x, 
the  total  normalized  potential  energy  of  the  system  is 


V„  • 4x2  - Vx, 


'1  2 


-2b- 


its  internal  energy  is  jx2,  and  the  Legendre  transform  of  this  is  the 
complementary  energy  having  the  value  ^b?2  in  equilibrium. 

Now  stepping  up  the  ladder  from  the  fold-type  potential  we  can 

regard  its  'partial  form'  -^x3  ♦ as  the  potential  energy  of  the 

perfect  version  of  the  stated  structure.  The  Legendre  transform  of  this 
is  the  complementary  energy  W of  the  perfect  version.  The  values  of 
this  when  the  imperfect  structure  is  in  equilibrium  have  the  same 
dependence 


(38) 

as  does  the  bifurcation  set  of  the  new  quartic  structure 
V4  * ^•x<*  ♦ 2t,2x3  * + 3Wx- 

This  fact  can  be  read  off  from  (26)  and  (35)  with  n = 3 and  X ■ 1. 

Geometrical  nonlinearities  frequently  induce  quartic  potentials  in  engin- 
eering structures,  and  if  we  change  the  origin  of  the  behaviour  variable 
here  by  writing  x = y - we  convert  this  potential  into 

v4  * -Jy4  - |v2[lto22  ' bl]  * ^^2 3 " ■§bib2  * 3W]  * constant. 

By  comparison  with  similar  expressions  given  by  Sewell  (1976)  for  compressed 
struts  and  shallow  arches,  we  can  identify  the  coefficients  here  with 
geometry  variables,  loads  and  imperfections  whose  buckling  values  are 
related  by  ( 38) . 
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By  constructing  the  Legendre  transform  of  the  internal  energy  (or  of 
the  ’perfect’  energy)  of  such  quartic  structures  with  linear  terms 
omitted,  we  can  arrive  at  the  bifurcation  set  of  another  structure  with 
swallowtail  potential,  and  so  on. 
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